Field Induced Staggered Magnetization and Magnetic Ordering in Cu 2 (C5Hi 2 N 2 )2Cl4 
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We present a 2 D NMR investigation of the gapped spin-1/2 compound Cu2(C5HioN2D2)2Cl4. Our 
measurements reveal the presence of a magnetic field induced transverse staggered magnetization 
(TSM) which persists well below and above the field-induced 3D long-range magnetically ordered 
(FIMO) phase. The symmetry of this TSM is different from that of the TSM induced by the order 
parameter of the FIMO phase. Its origin, field dependence and symmetry can be explained by an 
intra-dimer Dzyaloshinskii-Moriya interaction, as shown by DMRG calculations on a spin-1/2 ladder. 
This leads us to predict that the transition into the FIMO phase is not in the BEC universality 
class. 

PACS numbers: 75.10. Jm,75.40.Cx,75.300.Gw 



Since the pioneering work of Haldane for S=l anti- 
ferromagnetic ( AF) chains 0] , the existence of collective 
singlet ground states separated by an energy gap from 
the first triplet excited states has triggered a very large 
number of theoretical and experimental studies in low- 
dimensional quantum AF systems, as well as in 3D AF 
coupled dimers Q. Applying an external magnetic field 
H lowers the energy of the Ms = -1 component of the 
triplet band, inducing at some critical field value H c \ a 
quantum phase transition from a non-magnetic phase to 
a field- induced 3D long-range magnetic ordered (FIMO) 
ground state, which has been recently described in several 
cases as a Bose-Einstein condensation (BEC) of triplet 
excitations 0, S S 01- However, anisotropic terms in 
the spin-Hamiltonian, like Dzyaloshinskii-Moriya (DM) 
interactions or staggered g tensors, often open a gap 
at H c i, and change the universality class of the tran- 
sition 0, 111- Among numerous spin-liquid systems, 
Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 (Cu(Hp)Cl in short) has long been 
considered as the prototype of a S — 1/2 two-leg spin 
ladder in the strong coupling limit (Jj_ 3> J\\) Its 
phase diagram in the H-T plane consists of a gapped spin 
liquid phase below H c \ , a 3D ordered magnetic phase be- 
tween H c i = 7.5 T and H C 2 = 13 T, and a fully polarized, 
gapped phase above -ff c2 . In the intermediate phase, be- 
tween H c i and H c2 , the interactions between the "mag- 
netic ladders" lead to the FIMO phase [HE 03. closc 
to the first quantum critical point (H = H c \, T = 0), 
Tfimo(H) varies as (H -H c )^ a , where a was found to 
be close to 3/2, as expected in a BEC description W. 
However, inelastic neutron scattering measurements |l4j 
have shown that the exchange paths are not clearly iden- 
tified in this system and are more complicated than in 
the ideal picture described above. 

In this paper, we present a 2 D NMR study of a 



Cu(Hp)Cl single crystal (0.35x0.3x0.15 mm 3 size), in 
which the N-H groups have been replaced by N-D groups. 
2 D spectra have been recorded in the field range 3.5-15 T, 
and in the temperature range 50 mK-50 K. From the 
study of the quadrupolar couplings as a function of T 
and H, we rule out previous interpretations of the origin 
of FIMO as being magnetoelastic [l2l . We demon- 
strate the existence at low T of a sizable field-induced 
transverse staggered magnetization (TSM) which starts 
to grow several tesla below H c \ , remains nearly constant 
within the FIMO phase, and decays within a few tesla 
above i? c2 . Within the FIMO, it coexists with another 
TSM of different symmetry corresponding to the order 
parameter of the FIMO. Such a behavior has never been 
reported so far. We show that the TSM outside the 
FIMO is well explained by Density Matrix Renormal- 
ization Group (DMRG) calculations on a S = 1/2 spin 
ladder in the strong coupling limit if a DM interaction is 
introduced on the rungs. 

Selective deuteration provides several advantages over 
proton NMR. First, it divides by 6 the number of sites 
contributing to the spectra. There are only 8 inequiv- 
alent sites for an arbitrary orientation of H, and 4 for 
H || b or in the (a, c) plane, due to the symmetry of 
the space group P2i/c. Secondly, the substituted 1 H are 
the closest to Cu atoms. Half of them are involved in 
the presumed exchange path J«. These 1 H sites show 
an anomalous shift at low T, which has been previously 
considered as a hint for the mag netoelastic character of 
the transition into the FIMO [l3j. Finally, 2 D has a spin 
1 and a quadrupolar moment. Therefore, each inequiv- 
alent 2 D site gives rise to two lines. Their average posi- 
tion depends on the magnetic hyperfine shift K due to 
the coupling with the electronic spins borne by the Cu 2+ 
atoms through the hyperfine tensors. Their splitting Svq, 
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FIG. 1: (color online) T dependence of 2 D spectra for H = 
7.75 T || [111]. The quadrupolar doublet corresponding to 
the line A at the right hand side of the spectra is identified 
below 1.2 K. Its position is emphasized by solid lines, which 
are just a guide for the eyes. Below 600 inK, each component 
of the doublet splits into two lines (dashed versus solid line) , 
due to the entrance into the FIMO phase. Inset: Structure 
of Cu2(C 5 Hi2N 2 )2Cl4. The 10 protons attached to the C 
atoms are not shown. The other ladders can be deduced by 
translation and rotation around the 6 axis. 
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FIG. 2: a) T dependence of Suq for the site corresponding to 
the line A for H = 7.75 T and || to [111] b) H dependence 
of 5vq of the 4 sites at 50 mK for H | c*. In both cases, 
5vq remains constant. These data demonstrate the absence 
of any structural transition at the quantum critical fields H c i 
and H c2 at T ~ 0. 



due to the quadrupolar coupling with the environment, 
is very sensitive to atomic displacements or structural 
distortions. Both K and Svq strongly depend on the ori- 
entation of H with respect to the crystalline axes. So, 
from the simultaneous determination of 5vq and K ver- 
sus H and T, one can decide whether the variations of K 
are due to atomic displacements or not. As a drawback, 
2 D NMR is much less sensitive than 1 H NMR. All spectra 
were obtained by sweeping the frequency at fixed H value 
and summing the Fourier transforms of the echo [lij . As 
long as the hyperfine and the quadrupolar couplings are 
small perturbations with respect to the Zccman energy, 

5v l ± = 2 7 Shl ± 5^/2 (1) 

where 5h\ = K l H = ^ Z H. K* is the 

hyperfine shift of the deuteron site i, A 1 ' 1 the hyperfine 
field tensor between copper site I and deuteron i, while 
g' and \ l are respectively the g and the susceptibility 
tensors of copper I. Both A*' z and g' arc T independent 
in the absence of structural change. 

Before considering the T and H dependence of x(T), 
we first demonstrate that the quadrupolar couplings are 
T independent below 10 K and H independent at T ~ 
(50 mK). Fig. 1 shows 2 D spectra as a function of T for 
H = 7.75 T || [111]. In this particular orientation, the line 
A at the right hand side of the spectra remains isolated 
in the whole T range 50 K-50 mK. Several approaches 
have been used to determine the T dependence of 5vq 
for this site. Between 50 to 10 K, where all the 1 H lines 



have a similar T dependence, we made use of equation 
(1) by plotting 8va(T) versus : H shift 5^_f/(T), which re- 
flects directly the variation of xs{T). The extrapolation 
to 5v H = gives 5v% = 170 ± 8 kHz. Below 1 K, the 
second line of the doublet involving line A also becomes 
isolated in the spectra and the value of 8vq can then be 
directly determined. It is found constant through the 3D 
transition and its value 175 ± 6 kHz (Fig. 1 and 2a) is 
consistent with that determined by the high T procedure. 
These results strongly support the absence of structural 
transition in the whole temperature range 50 K-50 mK, 
and in particular at Tfimo- Let us now consider the H 
dependence of 5vq at low T. For that purpose, we used 
the spectra recorded at 50 mK in the field range 6-14 T, 
with H || c*. In this case the symmetry reduces the num- 
ber of inequivalent 2 D sites to 4, and it is possible to fit 
the whole spectra. Fig. 2b shows the values of 8vq for the 
4 sites as a function of H . They remain constant in the 
whole H range, in particular throughout the two quan- 
tum transitions at H c i=7.5 T and H C 2 ~ 13 T. This rules 
out the existence of a structural transition at T = 50 mK 
associated with the quantum magnetic phase transition. 
This corrects our previous interpretation on the nature 
of the FIMO 01 j an d contradicts the recently published 
X-ray experiments 0|. Our experiments also show the 
absence of any hysteresis upon entering the FIMO by 
varying H at fixed T = 50 mK. 

We can now focus on the T and H dependence of the 
shift K which directly reflects that of \- I n Fig. 3a, we 
have reported the T dependence of the hyperfine shift 
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3: a) Comparison between the T-dependence of A zz \ z 



FIG. 3: aj Uompa: 

(dashed line) and the shift K (•) of the line A at H 
7.75 T || [111], A zz = -470 G//i B . Below Tfimo = 600 mK, the 
line splits into two components, and solid squares correspond 
to their average position, b) The difference A zz Xzz(T) — K(T) 
(♦) is due to the onset of a transverse magnetization (see the 
text). The splitting S(T) (o) of the A line below Tfimo is 
proportional to the order parameter of the FIMO. 



K = Sh z /H of the A line and of the product of the lon- 
gitudinal spin susceptibility Yzz (T) by the hyperfine field 
A zz in the orientation [111] [lj|. These two quantities are 
no longer proportional below 3 K, and this can only be 
explained by the growing of a transverse magnetization: 



6h z (T) = A zz M z (T)+SH tr (T), 



(2) 



where SH tr = A ZX M X (T) + A zy M y (T). In Fig. 3b, the T 
dependence of the contribution of this transverse magne- 
tization to K(T) is plotted. Below Tfimo we consider the 
mean position of the split line (solid squares) and com- 
pare it to that of the splitting S which is proportional to 
the order parameter of the FIMO. 

More information on SH tr can be obtained by looking 
at its H dependence between 6 and 14 T at 50 mK, in the 
orientation H \\ c* . The experimental values of the shifts 
as a function of H arc shown in Fig. 4 and compared 
with their expected variations in the absence of 5H tr , i.e. 
A l zz M z {H). The solid and the dotted lines correspond to 
the extremal values of A l zz of the 4 hyperfine fields deter- 
mined between 50 and 10 K. One clearly sees before H c i 
the growing of two by two opposite values of SH tr , which 
remain nearly constant within the FIMO, and then start 
to decrease to zero. This immediately points to the ex- 
istence of a staggered magnetization, as in SrCu2(B03)2 
[2f1 |. However, it must be noticed that for both orienta- 
tions of H with respect to the crystal, the field induced 
TSM docs not split the NMR lines, as clearly seen outside 
the FIMO. This means that the number of inequivalent 
sites remains the same outside the FIMO (4 or 8 when 
H || c* or [111]). This differs from the effect of the spon- 
taneous TSM arising in the FIMO, which clearly splits 
them. The absence of splitting due to the field induced 
TSM implies that this TSM respects the full symmetry 
of the crystal space group, while that corresponding to 
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FIG. 4: a) H dependence of the hyperfine experimental shift 
of the 4 sites (▼, V, A) for H | c*. Between H c i and H C 2, 
the mean position of the split lines is plotted. The solid and 
the dotted lines are defined in the text, b) H dependence of 
the fraction of the shift due to the field induced TSM. 



the order parameter of the FIMO does not, as expected 
in a phase transition. 

The progressive appearance of a TSM upon approach- 
ing H c i at low T must come from interactions that break 
the SU(2) symmetry and mix the S = ground state with 
the lowest triplet S = 1 excitations. This can be due to 
DM interactions and/or staggered <?-tensors on dimers 
without inversion symmetry [20j , which is the case for 
the shortest Cui-Cu2 pairs |9|. To model the effects of 
DM interactions, we have performed DMRG calculations 
on a spin- 1/2 ladder assuming a staggered distribution of 
DM vectors T> on the rungs, as the simplest model con- 
sistent with the presence of an inversion center between 
neighboring dimers |2lj . Since the exchange paths arc 
still debated, other geometries including e.g. diagonal 
inter-rung couplings could in principle be realized. How- 
ever, preliminary results show that, as far as the TSM 
perpendicular to T> is concerned, the field dependence is 
essentially the same up to an overall factor |22j. As seen 
in Fig. 5a, the results reproduce the experimentally ob- 
served TSM fairly well for reasonable ratios T>/J± =0.05 
and J||/Jj_=0.2, with J± = 13 K. Furthermore, since m*j_ 
oc H x T> , it is easy to show, considering the symmetry 
transformations of a pseudo- vector in the group P2i/c, 
that the field induced TSM due to the DM interaction 
cannot split the 2 D lines, whatever the orientation of H: 
the alternation of the DM vector from one rung to the 
adjacent one corresponds to the symmetry of the lattice 
so that, if Si and 5*2 correspond respectively to Cui and 
Cu2 electronic spins, T> ■ S\ x S2 remains invariant. 

The situation is different for the TSM corresponding 
to the order parameter of the FIMO, which splits all the 
2 D lines whatever the field orientation. This means that 
it breaks the inversion symmetry between two adjacent 
rungs. An example of a possible pattern is shown in 
Fig. 5b. This would be consistent with a BEC in weakly 
coupled dimers systems into a uniform phase in which 
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FIG. 5: (color online) a) H dependence of the field induced 
TSM (V) and the fit (solid line) by the DMRG calculation 
with J ± = 13 K, J\\/J± = 0.2 and V/J ± = 0.05. The order 
parameter of the FIMO is proportional to S(H) (solid cir- 
cles), b) Schematic representation of the two "ladders" and 
the Cu(l) and Cu(2) sites viewed from the c* direction. Note 
an inversion symmetry between two successive rungs. Left: 
Schematic representation of the field induced TSM. Its direc- 
tion on one rung is fixed by H x T> and the whole pattern 
can be deduced by applying the symmetry operations of the 
lattice. Right: TSM pattern consistent with the order param- 
eter of the FIMO. The direction in the plane is unknown and 
fixed by the inter-dimer interactions. 

the TSM is the same on all rungs. However, the BEC 
description is here inadequate since, due to the field in- 
duced TSM, the isotropy in the plane perpendicular to 
H is lost. So, the ordered phase can only break a dis- 
crete symmetry, and the phase transition into the FIMO 
cannot fall into the universality class of BEC. This is 
actually supported by the present data. Indeed, T C (H) 
varies as (H — H c i) - 55 , consistent with the results re- 
ported in [r3|, but in slight disagreement with the BEC 
prediction (H - ff c i) 2/3 H3- Moreover, the field de- 
pendence of the order parameter within the FIMO, as 
revealed by the splitt ing S, was found to be (H — H c \)° 
with = 0.34 ± 0.04 [ljj, instead of the mean-field value 
(3 = 0.5 expected for a BEC in 3+1 dimensions. Un- 
fortunately, the exponent for the transition as a function 
of T at fixed H could not be reliably extracted from the 
present data. So, the precise determination of the nature 
of the phase transition into the FIMO is left for future 
investigation. Finally, the presence of DM interactions is 
normally accompanied by a gap opening, corresponding 
to some kind of level anticrossing |2jj. In Cu(Hp)Cl, a 
fast decrease of the nuclear spin-lattice relaxation T-j -1 
has indeed been observed below 1 K in the field range 
of the FIMO, even at H = 7.75 T for which T F imo = 
300 mK These data are consistent with the gap 



estimated by DMRG 

In conclusion, we have shown that in Cu(Hp)Cl two 
types of transverse staggered magnetization appear. The 
first one, which extends outside the 3D FIMO phase, is 
field induced and breaks the SU(2) symmetry, but not 
the crystal symmetry. We propose that it is due to the 
presence of a DM interaction T> on the strong dimers, 
and we have shown by an explicit calculation in a lad- 
der geometry that a quantitative fit can be achieved for 
a reasonable value of the parameters (T>/J± ~ 0.05). 
The second one, which corresponds to the order param- 
eter of the FIMO, breaks the inversion symmetry of the 
structure. We think that these features are generic for 
an assembly of interacting non centro-symmetric dimers, 
and we suggest that they will change the universality 
class of the field induced phase transition. 
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